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The most general Two Higgs Doublet Model potential without explicit CP violation depends on 
10 real independent parameters. There are two different ways of restricting this potential to 7 inde- 
pendent parameters. This gives rise to two different potentials, Via) and V(b)- The phenomenology 
of the two models is different, because some trilinear and quartic Higgs couplings are different. As an 
illustration, we calculate the decay width of h° — > 77, where precisely due to the different trilinear 
couplings the loop of the charged Higgs gives different contributions. We also discuss the possibility 
for the existence of a light fermiophobic Higgs. 
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I. INTRODUCTION 



Despite the great success of the standard SU(2) x U(l) electroweak model (SM), one of its fundamental principles, 
the spontaneous symmetry breaking mechanism, still awaits experimental confirmation. This mechanism, in its 
minimal version, requires the introduction of a single doublet of scalar complex fields and gives rise to the existence 
of a neutral particle with mass ran. The combined analysis of all electroweak data as a function of tub favors 
a value of run close to 100 GeV/c 2 and predicts with 95% confidence level an upper bound of mji < 200 GeV/c 2 . 
Hence, one can still envisage the possibility of a Higgs discovery in the closing stages of the LEP operation. 

Nevertheless, even if this turned out to be true, one still would like to know if there is just one family of Higgs fields 
or, on the contrary, if nature has decided to replicate itself. In our view this is the main motivation to consider multi 
Higgs models. In this paper we continue the study of the two-Higgs-doublet model (2HDM). Following our previous 
work ||, we examine models without explicit CP violation and which are also naturally protected from developing 
a spontaneous CP breaking minimum. There are two different ways of achieving this. To illustrate the different 
phenomenology we calculate, in both models, the decay width for the process h — > 77, which can be particularly 
relevant if h is a fermiophobic Higgs. 



II. THE POTENTIALS 



The Higgs mechanism in its minimal version (one scalar doublet) introduces in the theory an arbitrary parameter — 
the Higgs boson mass m# . In fact, the potential depends on two parameters, which are the coefficients of the quadratic 
and quartic terms. However the perturbative version of the theory replaces them by the vacuum expectation value 
v = 247 GeV and by vein . If we generalize the theory introducing a second doublet of complex fields, the number of 
free parameters in the potential V grows from two to fourteen. At the same time, the number of scalar particles grows 
from one to four. In this general form the potential contains genuine new interaction vertices which are independent 
of the vacuum expectation values and of the mass matrix of the Higgs bosons. However, these new interactions can be 
avoided if one imposes the restriction that V is invariant under charge conjugation C. In fact, if <E>i with i = 1, 2 denote 
two complex scalar doublets with hyper-charge 1, under C the fields transform themselves as — > cxp(iai)<I>* where 
the parameters a.i are arbitrary. Then, choosing a\ = oli = 0, and defining x\ = 4>\4>\, X2 = <f>\4>2, £3 = ^H{(j>[(j>2} and 
X4 = ^s{cf>\(f)2} it is easy to see that the most general 2HDM potential without explicit C violation^], is: 

V = -fJ-lxi - p? 2 X2 - H 2 2 X 3 + ^l x l + ^2x1 + A3X3 + A4IE4 + \5X1X2 + \qXiX 3 + X7X2X3 . (1) 

In general, the minimum of this potential is of the form 

«•*>-£(!.) <w 
<**>=7fO) ■ < 2b » 

in other words it breaks CP spontaneously. To use this potential in perturbative electroweak calculations the 
physical parameters that should replace the A's and /x's are the following: 

i) the position of the minimum, v±, V2 and 9, or alternatively, v 2 = v'f + v 2 , tan/3 = ^ and 9; 

ii) the masses of the charged boson m + and of the three neutral bosons mi, m,2 and m^; 

iii) and the three Cabibbo like angles a%, 0*2 and a.3 that represent the orthogonal transformation that diagonalizes 
the 3x3 mass matrixp] of the neutral sector. 



x At this level C conservation is equivalent to CP conservation since all fields are scalars. 

2 The mass matrix corresponding to the neutral components [T3 = — |J of the doublets is a 4 x 4 matrix, but one eigenvalue 
is zero because it corresponds to the Z would be Goldstone boson. 
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In a previous paper M we have examined the different types of extrema for potential V. In particular it was 
shown in |^] that there are two ways of naturally imposing that a minimum with CP violation never occurs. This, 
in turn, leads to two different 7-parameter potentials. The first one, denoted Via)-, is the potential discussed in the 
review article of M. Sher 3| and corresponds to setting (j,f 2 — ^6 = A 7 = in equation (|l|). The second 7-parameter 
potential, that we shall call Vib) > is essentially the version analyzed in the Higgs Hunters Guide ||| and it corresponds 
to the conditions \q = A 7 — and A3 = A4. As we have already pointed out M but would like to stress again, these 



potentials have different phenomenology. This is illustrated in section III when we consider the fermiophobic limit of 
both models. 

Since V(A) an d V(b) do not have spontaneous CP- violation, the number of so-called "physical parameters" is 
immediately reduced to seven. In fact, 9 — and only one rotation angle, a, is needed to diagonalize the 2x2 mass 
matrix of the CP-even neutral scalars. This is clearly seen if we transform the initial doublets <&i into two new ones 
Hi given by 



H 2 J ^Tvl V »i J\ $2 



(3) 



In this Higgs bases, only Hi acquires a vacuum expectation value. Then, the T3 — +A component and the imaginary 
part of the T3 = — | component of Hi are the and Z would be Goldstone bosons, respectively. The C-odd neutral 
boson, A , is the imaginary part of the T3 = — J component of H 2 . On the other hand, the light and heavy CP-even 
neutral Higgs, h° and H°, are linear combinations of the real parts of the T3 = — | component of Hi and H 2 . 

Notice that V(a) is invariant under the Z 2 transformation $i — > $! and $2 — > —&2, whereas in V^b) onr y the 1A.2 
term breaks the ?7(1) symmetry, $ 2 e la &2- Because this breaking occurs in a quadratic term it does not spoil the 
renormalizability of the model. Hence, in both cases the terms that were set explicitly to zero, will not be needed 
to absorb infinities that occur at higher orders. The complete renormalization program of the model based on V^a) 
was carried out in ||. The results for Vm) are similar but the cubic and quartic scalar vertices have to be changed 
appropriately. 

For the sake of completeness we will close this section with a summary of the results that will be used later. As we 
have already said they are not new and can be obtained either from || or j| . We agree with both. 
For V(a) the minimum conditions are 

= T 1 =v 1 (-^ 2 1 + X 1 v 1 2 + X + v 2 2 ) (4a) 
= T 2 = v 2 ( -i4 + A 2 «2 2 + A+ Vl 2 ) (4b) 

with A+ = I (A3 + A5). They lead to the following solutions: 
either i) 



A 2 ^i — A + /i 2 
A1A2 — A^_ 

A1M2 ~ 
A1A2 — A + 



(5a) 
(5b) 



ii) 



v\ = (6a) 
, 2 2 = # . (6b) 



The masses of the Higgs bosons and the angle a are given by the following relations: 



■in 



2 

H+ 



m 2 H o 



-A 3 (v 2 + v 2 ) (7a) 
%> = i(A 4 -A 3 ) {vl + vl) (7b) 

h o = Xiv 2 + \ 2 v\ ± yj (\ lV 2 - \ 2 v 2 ) 2 + v 2 v 2 (\ 3 + A 5 ) 2 (7c) 
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tan 2a = 



v 2 vi (A 3 + A 5 ) 



Xivf - X 2 v 2 

On the other hand, for V^b) t ne minimum conditions are 

n rp lA.2 

= Ti — - v 2 

= T 2 — v i 

with the Ti given by the previous equations (4). The solution of this set of equations is 



(8) 

(9a) 
(9b) 



Ai - A 2 ± sj (Ai - A 2 ) 2 - 4 (Ai - A+) (A 2 - A+) [(A+i, 2 - (A 2 t; 2 - ^) - ±^ 



2 (Ai - A+) (A 2 - A+) 



A 2 - Ai ± yj (Ai - A 2 ) 2 - 4 (A 2 - A+) (X 1 - A+) [(A+i; 2 - M |) (Ai^ 2 - M?) 



4^12 



2 (Ai - A+) (A 2 - A+) 



(10a) 



(10b) 



Notice that, in this case, the solution with vanishing vacuum expectation value in one of the doublets is not possible. 
Now the masses and the value of a are given by 



-A3 (vf + v 2 2 ) + ii\ 



2 u l 
2 



v? + vl 



viv 2 



2 _ 1 2 "l + ^ 



l 2 ,o j ,o=A lU 2 + A 2 i; 2 2 + i M 2 2 (^ + ^) 

± y (a iW 2 - A 2 « 2 + ± M 2 2 - ^)) 2 + (z^ 2 (A 3 + A 5 ) - ±^ 2 ) S 



(11a) 

(lib) 
flic) 



tan 2a 



2^iu 2 A + 



2^12 



Aiu 2 - A 2 wf 



±/v 2 

4P12 



I V2_ _ VI 
\ VI V2 



(12) 



III. THE FERMIOPHOBIC LIMIT 



Despite the fact that Vm) and V(b) are different, it is obvious that the gauge bosons and the fermions couplings to 
the scalars are the same for both models. In particular, the introduction of the Yukawa couplings without tree-level 
flavor changing neutral currents is easily done extending the Z2 symmetry to the fermions. This leads to two different 
ways of coupling the quarks and two different ways of introducing the leptons, giving a total of four different models, 
usually denoted as model I, II, III and IV (cf. e.g. ||). 

In here, we use model I, where only $ 2 couples to the fermions. Then, the coupling of the lightest scalar Higgs, h°, 
to a fermion pair (quark or lepton) is proportional to cos a. As a approaches 5 this coupling tends to zero and in the 
limit it vanishes, giving rise to a fermiophobic Higgs. 

Examining equations (^|) and (Jl2] ) we see that the fermiophobic limit (a = ^) can be obtained in potential A in 
two ways: either A+_= or V\ = 0. In potential B there is only one possibility 2uiw 2 A + = ^/i 2 2 . In this latter case, 

(13a) 

m 2 A0 + 2(A 2 - A+) v 2 sin 2 (3 (13b) 
m\ Q - 2(A+ - Ai)v 2 cos 2 /3 . (13c) 



equations (11) and ( |l2| ) give immediately: 

m\ n = 2 A+ (vj + v 2 2 ) 
m 2 H o = 2X 2 v 2 + 2X+v 2 = 
m^o = 2 Aid 2 + 2 X+v 2 , 

In the former case (V^), A + = gives 
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m 2 H0 = 2\ 2 v\ (14a) 
mlo=2\ 1 v 2 1 (14b) 

while vi = gives a massless h°. In this analysis we have assumed that V\ < i>2- The reversed situation leads to 
similar conclusions since one is then interchanging the role of the two doublets. 

The triple couplings involving two gauge bosons and a scalar particle like, for instance Z^Z^h , are always propor- 
tional to the angle S = a — j3. In particular, the couplings for h are proportional to sin 8 whereas the corresponding H° 
couplings are proportional to cos <5. This general results can be understood if one recalls the argument about the role 
played by the neutral scalars in restoring the unitarity in the scattering of longitudinal VF's, i.e. in 

w l W l -> w l W l- 

The restoration of unitarity requires that the sum of the squares of the W + W~h° and W + W~ H couplings adds up 
to a constant proportional to the SU(2) gauge coupling, g. 

Current searches of the SM Higgs boson at LEP put the mass limit at 89 GeV/c 2 [§. Since the production 
mechanism is the reaction e + e~ — > Z* — > Zh°, this limit can be substantially lower in the 2HDM if sin<5 is small. In 
our numerical application to the two 7 decay of a light fermiophobic h° we will explore the region sin 2 8 < 0.1 0. 

Bounds on the Higgs masses have been derived by several authors ||. Recently next-to-leading order calculations 
H in the SM give a prediction for the branching ratio Br(B — ► X s ~f) which is slightly larger than the experimental 
CLEO measurement |l(J. In model II the charged Higgs loops always increase the SM value. Hence, this process 
provides good lower bounds on m H ± as a function of tan (3 j^] . On the contrary, in model I the contribution from the 
charged Higgs reduces the theoretical prediction and so brings it to a value closer to the experimental result. This 
reduction is larger for small tan/3, since in model I the H + coupling to quarks is proportional to tan -1 (i. However, a 
small tan/3 gives a large top Yukawa coupling which leads to large new contributions to Rb, the Bq — Bq mixing. A 
recent analysis by Ciuchini et al. || derives the bounds tan/3 > 1.8, 1.4 and 1.0 for m H ± = 85, 200 and 425 GeV/c 2 , 
respectively. 

The Higgs contribution to the /9-parameter is PJ : 

Ap = 1 [sin 2 6 F(m 2 H ± , m 2 A0 , m% ) + cos 2 6 F(m 2 H ± , m 2 A0 , mf o)] (15) 
Vq-k v 

where 

. be b ab a ac a 
F(a, b. c) = a + In - — In — — In — 

— c c a — b b a— c c 

Since the current experimental value of p = 1.0012 ± 0.0013 ± 0.0018 (l^] exceeds the SM prediction by 3er, one 
should at least try to avoid a positive Ap.^|. A simpler examination of the function F(a, 6, c) shows that this is 
impossible if m H ± is the largest mass. On the other hand, if m^o > m H ± one obtains a negative value for Ap 
which grows with the splitting m A o — m H ±. In line with our limit (sin 2 <5 < 0.1), negative values of Ap of the 
order of the experimental statistical error, i.e. Ap w — 10~ 3 , can be obtained essentially in two ways. Either with 
a large m#± ~ 300 GeV/c 2 but with a modest m A o — m H ± splitting (m A o w 340 GeV/c 2 ) or with a smaller 
m H ± w 100 GeV/c 2 but with 771,40 w 200 GeV/c 2 . The variation of Ap with m} x a is rather modest, less than 10% for 
the range 20 GeV/c 2 < m h o < 100 GeV/c 2 . With seven parameters in the Higgs sector it is difficult and not very 
illuminating to discuss in detail all possibilities. So, this discussion should be regarded as a simple justification for 
the fact that a fermiophobic Higgs scenario is not ruled out by the existing experiments. We would like to stress, 
that there could exist a light h° almost decoupled from the fermions (asi I) and at the same time with a small LEP 
production rate via the Z-bremsstrahlungs reaction Z* — > Zh° (sin 2 S rs 10 _1 ). If such a boson exists it will decay 
mainly via the process h — > 77. 



IV. THE DECAY h° 



77 



The decay h° — > 77 is particularly suitable to illustrate the fact that Via) an d Vm) gi ye rise 1° different phe- 
nomenologies. In fact, the decay occurs at one-loop level and for a fermiophobic Higgs one has vector bosons and 
charged Higgs contributions. The latter are different for models A and B, because the h°H + H~ vertex is different. 
It is interesting to point out how this difference arises. Since the term in A4 does not contribute to this vertex, both 
potentials give rise to the same effective h°H + H~ coupling, gh°H+H- 1 namely: 



a A more recent SM fit gives p = 0.9996 + 0.0031(-0.0013) |§|. 
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[h°H + H ] = 2 V2A2 cos 2 (3 cos a + V2X3 sin a cos (3 sin (3 — v\ A5 cos 2 /3sina 
—2 W1A1 sin 2 (3s'ma + V2A5 sin 2 /3 cos a — U1A3 cos a cos /3 sin /3 



(16) 



However, as we have already said, what is relevant for perturbative calculations is the position of the minimum of 
V and the values of its derivatives at that point. This means that one has to express all coupling constants in terms 
of the particle masses. This is simply done by inverting equations (7) and (11). The result is 



[h°H+H-i- - --^ 



cos (a+(3) 



(A) ~^l m ' l °" sin 2(3 



— (m 2 H+ — \m 2 ha ) sin (a — (3) ^ 



and 



[h°H + H-] [B) = 



fJ 



2 2 \ cos(a+/3) 
m h0 - m A0 ) 



mw \ v '" " ' sin 2(3 
which clearly shows the difference that we have pointed out. 



(m 2 H+ - inijo) sin [a- pi) 



(17) 



(18) 





(q 10-q20,-q21) (q10,0) 
(q20,q21) > h 




(q 10-q20,-q21) (q10,0) ^ 

(q20,q21l > 




{q 10-q20,-q2 1 ) (q10,0) 
(q20,q21) > 




{q 10-q20,-q2 1 ) (q10,0) 
{q20,q21) > 





(q 10-q20,-q21) (q10,Q) 
(q2Q,q21l > 




{q 10-q20,-q2 1 ) (q10,0) 
{q20,q21) > 




{q 10-q2O,-q2 1 ) 
(q20,q21l > 




{q 10-q20,-q2 1 ) 
{q20,q21) > 




(q1D-q20,-q21) (q1D,0) 



(q10-q20,-q21) (q10,0) 



(q10-q2O,-q21) (q10,0) 







(q10-q20,-q21) (q10,0) 



{q10-q20,-q21) (ql^O) 



{q10-q20,-q21) (q10,0) 



G 



jqJOflJI^ 






(q 10-q20,-q^1) (q1D,0) 





{q 10-q2Q,-q2J) (q10,Q) 





{q 10-q2O,-q2J) (q10,Q) 
jq2£fl21^_ 



{q 10-q2Q,-q2J) 
(q2°,q21) > 




(q1D-q20,-q21) (q1D,0) 



(q10-q20,-q21) (q10,0) 



(q10-q2O,-q21) (q1D,0) 




FIG. 1. The contributing graphs to h° — > 77 in the fermiophobic limit. 



In Fig. [| we show all the diagrams that were included. A previous work by Diaz and Weiler ^| did not include 
the Higgs-bosons diagrams. Our calculation, in the 'tHooft-Feynman gauge, was done with xloops |l^Jl^]. We have 
been using this program to calculate other amplitudes in the framework of the 2HDM [jlTj . Throughout this process 
we have made several checks on the computer results. In this particular case we have verified that the contribution 
of the vector boson loops agrees with a calculation done by M. Spira et al. Jl8| using the supersymmetric version of 
the 2HDM. 

In Fig. ^| we show the product m h o times the decay width (r) for the process h° — > 77 in model A as a function of 
5 and for several values of m^o and a fixed value of m H + . This function shows a gentle rise with m^o which reflects 
the proportionality between gh°H+H- an d Tn\a- Looking at this coupling constant one could naively assume that 
there would be an enhancement for j3 approaching 7r/2, i. e., in our plot, when S approaches zero. However, a close 
examination shows that such an enhancement does not exist. On the contrary, the coupling vanishes in this limit, 
since m^o goes to zero when f3 — > tt/2. Alternatively, if one keeps m^o fixed, then the mass relation 



m h o 



(19) 



imposes a lower bound for (3. In Fig. || the dotted line gives this limit, evaluated assuming Ai ~ 1/2. The dashed 
area shows the exclusion region implied by the LEP experimental results. In the work of Ackerstaff et al. jl9| an 
experimental bound on the SM h — > 77 branching ratio is derived. For a fermiophobic Higgs with < mw the 77 
branching ratio is one. On the other hand, the production mechanism is suppressed by a factor sin 2 5. Hence, we 
have turned the OPAL experimental bounds into a bound on S. Fig. ^| gives the equivalent information for potential 
B. 
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Decay width of h 
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|5| 

FIG. 2. Dependence on 5 and m h o at 772^+ = 131 GeV and a = ^ in potential A. 



Decay width of h 




FIG. 3. Dependence on 8 and m h o at m H + = 131 GeV and a = ? in potential B. 

In Fig. |we plot, as a function of m h a, the ratio R, of the widthes calculated with potentials Vb and Va, respectively. 
According to the fermiophobic limit, we set a = ir/2 and 5 = 0.29. For the other relevant masses we have used 
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m H + — 200 GeV/c 2 and m A o = 250 GeV/c 2 . In the range of variation of m h a, i.e., 20 GeV/c 2 < m h o < 120 GeV/c 2 , 
R decreases smoothly from 25 till 3. However, it is misleading to assume that potential A always gives smaller results. 
This is clearly shown in Fig. |5| where we plot the same function R evaluated with the same parameters except for 
m^a that was set up to 120 GeV/c 2 . Again, R is a decreasing function of m^o that has a zero for m h a around 
70 GeV/c 2 and increases afterwards. However, in this case, the values obtained with potential B are smaller than the 
corresponding ones for potential A. 



Ratio of h° decay widthes 




I 1 1 1 1 1 1 

20 40 60 80 100 120 

m(h) [GeV] 

FIG. 4. Ratio of the decay widthes from V(b)/V(a) with 8 = 0.29 and m H + = 200 GeV and m A o = 250 GeV. 

This behavior can be qualitatively understood if one examines the coupling constants [h° H + H~\i^\ and 
[h°H + H~]( B ) given by equations fll7| ) and (]lSj), respectively. In the range of m h o that we are considering, and 
for the same values of a, and m H + , the coupling corresponding to potential A is always negative and decreases 
from about —85 GeV/c 2 till —230 GeV/c 2 . On the contrary, the coupling constant corresponding to potential 
B is positive. For large values of m^o (around 250 GeV/c 2 ), it decreases from 930 GeV/c 2 till 780 GeV/c 2 for 
20 GeV jc 2 < myfi < 120 GeV/c 2 . These values of the coupling constant, when compared with the corresponding ones 
for potential A, explain the qualitative behavior of the ratio R given in Fig. |J. The explanation of Fig. || is more 
subtle, but again, it depends on the coupling constant of potential B. In fact, when m^o — 120 GeV/c 2 the coupling 
corresponding to potential B starts at 100 GeV/c 2 and decreases smoothly till —60 GeV/c 2 , having a zero around 
to/jo = 95 GeV/c 2 . This behavior has two consequences. When the coupling is positive, its order of magnitude is the 
correct one to almost cancel the W-loops contributions to the width. Hence, R is small because potential B gives a 
small width. This cancellation is exact for m h o around 70 GeV jc 2 and after that, because the coupling changes sign, 
the charged Higgs contribution adds up to the normal W-loop result. Hence R increases. 
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Ratio of h decay widthes 




m(h) [GeV] 

FIG. 5. Ratio of the decay widthes from Vm)/V(A) with S — 0.29 and m H+ — 200 GeV and m A o = 120 GeV. 

Despite the fact that the [/iWW] coupling is suppresed by sin 5, one should keep in mind that when m,h is larger 
then mw the decay channel h —* WW* — > Wqq starts to compete with the 77 channel. We have evaluated the 
WW" decay width and in table | we show some results in comparison with the width for the 77 channel evaluated for 
potential A and mu+ = 100 GeV . The table is representative of a situation that can be summarized qualitatively as 
follows: i) for small S (S = 0.1) the WW* width is comparable with the 77 width for m/, = 120 GeV/c 2 ; ii) for large 
6 (6 = 0.3) even at m/j = 120 GeV/c 2 the WW* decay width is already larger than the 77 width by a factor of ten. 







S = 0.1 






S = 0.3 






h -» WW* 




h -> 77(A) 


h -> WW* 




h -» 77(A)' 


90 


2.6 x lO -13 




0.4 x lO -3 


2.3 x 10 _b 




0.6 x 1Q-' A 


120 


2.2 x 10~ 3 




2.7 x 10~ s 


1.9 x 10 -2 




2.2 x 10~ 3 


150 


5.7 x 10~ 2 




1.5 x 10" 2 


5.0 x 10" 1 




8.0 x 10" 3 



TABLE I. Comparism between the widthes for the WW* and 77 channels. 
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V. CONCLUSION 



We have examined the 2HDM where the potential does not explicitly break CP violation and furthermore it is 
naturally protected from the appearance of minima with CP violation pL There are two ways of accomplishing this, 
leading to two different potentials Va and Vb ■ Va is invariant under the discrete group Z2 and Vb is invariant under 
U(l) except for the presence of a soft breaking term. These two symmetries ensure that the parameters that, at 
tree-level, were set to zero, are not required to renormalize the models. 

The potential Va and Vb have different cubic and quartic scalar vertices. Then, it is obvious that they give different 
Higgs-Higgs interactions. However, even before one is able to test such interactions, one could still sense these two 
different phcnomcnologies via Higgs-loop contributions. 

To illustrate this point we have considered a fermiophobic neutral Higgs, decaying mainly into two photons. The 
widthes for the decays calculated with both potentials can differ by orders of magnitude for reasonable values of 
the parameters. Clearly, with four masses and two angles as free parameters, it is not worthwhile to perform a 
complete analysis. Nevertheless, we believe that the results presented here are sufficient for illustrative purposes. The 
experimental searches in this area should be made with an open mind for surprises. 
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